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The Lorentz Integral Transform approach allows microscopic calculations of electromagnetic reac¬ 
tion cross sections without explicit knowledge of final state wave functions. The necessary inversion 
of the transform has to be treated with great care, since it constitutes a so-called ill-posed problem. 
In this work new inversion techniques for the Lorentz Integral Transform are introduced. It is shown 
that they all contain a regularization scheme, which is necessary to overcome the ill-posed problem. 
In addition it is illustrated that the new techniques have a much broader range of application than 
the present standard inversion method of the Lorentz Integral Transform. 

PACS numbers: 25.30.-c Lepton induced reactions - 25.30.Fj Inelastic electron scattering to continuum - 
25.20.Dc Photon absorption and scattering - 02.30.Uu Integral Transforms 


I. INTRODUCTION 


About a decade ago the Lorentz Integral Transform (LIT) method has been proposed in order to perform ab initio 
calculations of electroweak reactions with nuclei into the continuum Q. The great advantage of the method lies 
in the fact that a calculation of continuum wave functions is not necessary. Indeed the LIT approach reduces a 
continuum state problem to a bound-state like problem and thus it is sufficient to use bound-state methods. Hence 
it is not surprising that the LIT method has been applied to microscopic calculations of quite a few electroweak cross 
sections of various nuclei ranmng from A=2-7: inclusive electron scatteriM (see e.g. Hi): total photoabsorption 
cross sections (see e.g. 0, i Q): exclusive electromagnetic reactions 013]: photomeson production 0, and weak 
processes |in). The list of application shows that the LIT approach has proven to be an important progress opening 
up the possibility to carry out microscopic calculations not only for reactions of classical few-body systems (deuteron, 
three-body nuclei), but also for reactions of more complex nuclei. 

The LIT method proceeds in two steps. First one calculates an integral transform of the searched function i?(w) 
with a Lorentzian-shape kernel: 


A(crR, cri) 



+ erf ’ 


(Tr , (Tl > 0 . 


( 1 ) 


Then, in a second step, one inverts the obtained integral transform. The inversion, however, is in principle an ill- 
posed problem and the resulting consequences were already discussed in m- In order to introduce the reader to 
the problem we summarize the main points of that discussion. The inversion of o is unstable with respect to high 
frequency oscillations H. Adding such a high-frequency term AqR(uj) to R{oj) leads to an additional AoL(ctr, cti) 
in the transform. With growing H on the one hand and a constant non negligible oscillation amplitude AqR on the 
other hand AqL becomes increasingly small and might become even smaller than the size of errors in the calculation. 
Thus AqR could not be discriminated. If one reduces the error in the calculation one can push the not discriminated 
AqR to higher and higher Cl. However, even if the excluded frequency range is physically unimportant, one cannot 
simply find a solution of the response by application of the inverse operator, since the unphysical oscillations cannot 
be separated from the solution. As further explained in El one then has to use a regularization scheme for the 
inversion in order to avoid such problems. Therefore a specific inversion method with a built-in regularization was 
recommended. In fact this method has been used in all the above mentioned LIT calculations leading to very safe 
inversion results. On the other hand in all these calculations the various LITs have (i) a rather simple structure, where 
essentially only a single peak of R{uj) has to be resolved or (ii) a more complicated structure, which however could 
be subdivided in a sum of simply structured responses, where the various LITs have been inverted separately. Case 
(ii) was already encountered in [y, where the inclusive longitudinal deuteron electron response R{u!) was calculated 
at a constant momentum transfer. There it was necessary to separate out the Coulomb monopole and quadrupole 
transitions, which lead to a shoulder of the corresponding i?(a;) at the break-up threshold, while the rest of the 
response shows the typical quasi-elastic peak structure. Obviously it would be advantageous to invert the total LIT 
by just one inversion with a method appropriated for more complicated structures. In addition it is not guaranteed 
that one can always make the above subdivision in a sum of simply structured responses, e.g. in case of two rather 
close peaks in the same multipole transition. In such cases the standard inversion method would not be sufficiently 
precise. 
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The aim of the present paper is twofold. The main purpose is the investigation of the problem of the inversion of a 
more structured R{uj) and to introduce alternative inversion methods. As already pointed out the LIT method reduces 
a continuum state problem to a bound-state like problem, which can be solved with typical bound-state techniques. 
In recent years many different solutions for the bound-state problem of nuclei with A > 3 have been developed as 
illustrated for example in a rather recent paper about the ^He ground-state solution with a realistic nucleon-nucleon 
potential, where seven different theoretical methods were applied ^2- Many of these methods could also be used 
for a LIT calculation and as a second aim we want to show to potential future practitioners of the LIT method that 
reliable inversion techniques exist even in case of more structured response functions. 

The paper is organized as follows. In section II we give a brief outline of the LIT approach. The standard inversion 
method and alternative approaches are described in section III. Results with the various inversion techniques are 
discussed in section VI. 


II. THE LORENTZ INTEGRAL TRANSFORM METHOD 

In this section we give a brief outline of the LIT technique discussing however the approach for inclusive reactions 
only. The calculation of exclusive reactions is somewhat more complicated, but one proceeds in principle in a very 
similar manner 0,0 . 

As pointed out in the introduction the starting point of the LIT method P| is the calculation of the integral 
transform of R{uj) given in eq. (^. The function R depends on the internal excitation energy uj = Ef — Eq of a given 
particle system and contains information about the transition of the system from the ground state |5'o), with energy 
Eq, to the final states |'I'/}, with energy Ef, induced by an external probe. In case of an inclusive reaction i?(w) 
denotes the response function 

i?H= J d'i’f\{'i>f\d\'i>o)f6{Ef-Eo-Lj), (2) 

where O is a transition operator which characterizes the specific process under consideration. The response function 
and in principle also the LIT are observable quantities. However, in experiment R{uj) can be determined more directly, 
while L(crR,(Ti) is only accessible via the explicit integration of an experimentally over a sufficiently large ta-range 
determined Thus also in a theoretical calculation it is preferable to eventually determine i?(w). The key point 

of the LIT method is nonetheless a direct calculation of L((TR,cri), i.e. without explicit knowledge of i?(w). Only in 
a second step R is obtained from the inversion of the transform. The great advantage of the method lies in the fact 
that the generally very complicated calculation of final state wave functions jT/) can be avoided as will be discussed 
below. On the contrary a conventional calculation of R{uj) can only be carried out with the explicit knowledge of 

Given an Hamiltonian H of the particle system under consideration one may use the completeness of the eigenstates 
of E[ to show that L((Tr,(Ti) is determined via the following differential equation P] 

=6|^o) (3) 

with 

(T = a's. + i(Ti. (4) 

Since H has a real eigenvalue spectrum the corresponding homogeneous differential equation has only the trivial 
solution. Thus the solution of © is unique. This solution leads directly to the transform, in fact one finds 

A((TR,cri) = (4'|4'). (5) 

Different from a Schrodinger equation at positive energies one has a very simple boundary condition for the solution 
of @. Due to the localized source at the right hand side of 0 and because a is complex, ^ vanishes at large 
distances similar to a bound-state wave function. Therefore one can apply similar techniques as for the calculation of 
a bound-state wave function. It is evident that the much more complicated continuum state problem is completely 
avoided by the LIT method. 


HI. INVERSION METHODS FOR THE LORENTZ INTEGRAL TRANSFORM 

In order to invert the LIT one has to calculate L((TR,cri) for a set of ctr, values at fixed tji > 0. If one wants to 
determine R(oj) in a range uJmin < a; < Wmax one should choose ctr in the interval Wmin — o’! < ctr < Wmax + For 
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such ctr values the kernel in o takes a resonant form with the width determined by cri. The principal difhculties of 
the LIT inversion were already discussed in the introduction, hence in the following subsections we directly describe 
the standard LIT inversion method (subsection A) and the new inversion approaches (subsections B-D). 


A. Inversion with a set of basis functions 


The present standard LIT inversion method consists in the following ansatz for the response function 

N 

R{^') = X! ^riXn{uj', ai) (6) 

n—1 

with Lo' = uj — wth, where Wth is the threshold energy for the break-up into the continuum. The x„ are given functions 
with nonlinear parameters ai. A basis set frequently used in the LIT calculations of [1-11] is 

X„(w,ai) = exp(- —). (7) 

nuj 

In addition also possible information on narrow levels, like e.g. elastic transitions, could be incorporated easily into 
the set Xrt- 

Substituting such an expansion into the right hand side of m one obtains 

N 

L(CTR,Cri) = ^ C„Xn(CTR,Cri,Q!i) , (8) 

n—1 


where 


Xn(o’Rj O’! 


= r 

JO 


dio' 


Xniuj',ai) 

{lo' - crR)2 -I- af 


(9) 


For given ai the linear parameters Cn are determined from a best fit of L((Tr, (Ti) of eq. (|H1) to the calculated L((Tr, cti) 
of eq. for a number of ctr points much larger than N. If one takes the basis set o one should vary the nonlinear 
parameter 02 in a rather large range, while one can determine ai from the possibly known threshold behaviour of R 
or one can vary ai within a reasonable range. 

For any values of N, ai, and a 2 the overall best fit is selected and then the procedure is repeated for N = N + 1 till 
a stability of the inverted response is obtained and taken as inversion result. A further increase of N will eventually 
reach a point, where the inversion becomes unstable leading typically to random oscillations. The reason is that 
L(tTR, CTi) of eq. ||SJ) is not determined precisely enough so that a randomly oscillating R{lo) leads to a better fit than 
the true response. It is evident that the number of functions N plays the role of a regularization parameter and 
has to be chosen within the above mentioned stability region. Normally such a stability region is reached without 
greater problems, however, it can in principle happen that one does not find a stable result. Then one can either try 
to improve the precision of the calculated L(crR, ui) or use different basis sets. In case that the response exhibits an 
unexpected structure it is useful to decrease the parameter ai in order to have a better resolution in the transform. 

We should mention that the completeness of the basis set Xn is not really relevant. In fact any basis set is only used 
up to a relatively low value of N. Information which can be parametrized only via Xn with n > N is lost anyway. 
Therefore it is important to work with a basis set, where a relatively small N leads to a high-quality result of R{oj). 
Because of the variation of the nonlinear parameters ai of Xn many different basis sets are used and it is normally no 
problem to find a proper basis set. 

As already pointed out the standard LIT inversion method leads to high-precision results in case of simply structured 
response functions, but can have problems in case of more complicated structures. If one considers however just the 
low-energy part of the response the standard method is very reliable. This is due to the fact that the break-up 
threshold Wth is directly incorporated in the inversion {R{lo) = 0 for w < Wth) and also a known low-energy behaviour 
can be directly implemented {e.g. by the parameter ai in 0 ). In addition the inversion can be restricted to a small 
energy interval above threshold avoiding problems arising from structures at higher energies (see also [ 13 ). 


B. Inversion with wavelets 

The inversion with wavelets is formally very similar to the inversion with a set of basis functions described in the 
preceding subsection. One starts from the same expansion 0, but one uses wavelets instead of the basis functions 
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TABLE I: M and N^. values used for the wavelet inversion of section IV 


M=Q 

No 15 31 63 127 255 

M=1 

No 7 15 31 63 127 255 15 31 63 127 255 

Ni 7 15 31 63 127 255 7 15 31 63 127 

M=2 

Wo 7 15 3l 63 127 31 63 127 255 

Ni 7 15 31 63 127 15 31 63 127 

N2 7 15 31 63 127 7 15 31 63 


Xn- Here we will take the Mexican hat wavelet 

V3(a:) = (l-x^) exp(-y), (10) 

which is shown in fig-H We expand R{uj) in the range Wth < w < Wmax as follows 

M Nm 

R(^Co) = ^ ^ ^ tr) (H) 

m—Qn—1 

with 


= (12) 

where a is a nonlinear parameter and the Nm grid points ujm,n are equally distributed (uJm,n = n(uJma.x—0Jth) / (Wn+l))- 
Such an expansion leads to the following LIT 


M Nm 

L((^R, cri) = EE Cm,nV^m,n (^R; tt) 

m—0 n—1 


with 


^m,n(rrR, <Jl, o) 





(W - (Tr)2 + 0-2 


(13) 


(14) 


The wavelets ipm,n are strongly localized functions about oj = U!m,n in contrast to the basis functions Xn in subsection 
A, which generally are considerably different from zero over a rather large oj-range. The widths of the wavelets can 
be modified with the parameter a. In addition with the parameter M one can choose the presence of various scales 
of wavelets. Because of these wavelet characteristics one easily understands that they are much more appropriated to 
represent complicated structure of R{oj) than the basis functions Xn of the previous subsection. 

For the actual inversion problems discussed in the next section we use up to three different wavelet scales and 
various values for the Nm sets (see table I), in addition we include an additional basis function ipo, namely a constant 
((/3o = Co). For any of the 25 Nm sets of table I we proceed as follows. The nonlinear parameter a is varied over a 
wide range and for any value of a a best fit is performed leading to the determination of the linear parameters Cm,n- 
For any Nm set the overall best fit is then chosen. Of these 25 inversion results we choose again the overall best fit. 
In case that unrealistic oscillations are present in the inversion the next best inversion result is taken. Of course such 
a procedure is nothing else than a regularization scheme. 


C. Fridman approach applied to the LIT inversion 

The following method to invert the Lorentz integral transformation is based on the approach of V. M. Fridman 
la This method solves by iteration the Fredholm integral equation of first order 


dyK{x,y)f{y), 


F{x) = 


(15) 
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which has a solution and converges in the Hilbert space of square integrable real functions on [a,b] {L 2 {a,b)) if the 
kernel K is strictly positive definite and symmetric. The iterative equation is given by 


fi+iix) = fi{x) + X 


F{x) - f dyK{x,y)fi{y) 

J a 


^eNn 


(16) 


with /o G ^ 2 ( 0 , b) and 0 < A < 2Ai, where Ai is the smallest characteristic number of K. Using the spectral theorem 
it is evident that the iteration of eq. cni) diverges if at least one eigenvalue of if is in ] — 00 ,0]. 

To apply this method to the inversion of the Lorentz integral transform 


POO 

L{(Jb.,cfi)= / dujK{a^,ai,uj)R{uj) , if (ctr, cti, w) = 


(W - (Tr) 2 + CTj^ 


(Tl > 0 


(17) 


one has to discretisize the integration and symmetrize the integral kernel if. In practice the integration in eq. 03 
ranges from Wth to Wmax + ^i, since contributions beyond Wmax + ci are of minor importance. The numerical integration 
is done using the simple formula: 


dxf{x) 

h 


Xi 


n 

'^Wif{Xi), 

i=0 

b-a h , ■ n n 

-, Wo = Wn = - , Wi = h, I e |1,... ,n - 1} , 

n 2 

a + hi, i G {0,...,n} 


(18) 


Therefore the symmetrized iterative equation is given by: 


ii;+i(a;i) = i?/(a;i) + A 


ifi"‘^ = 


WiWj 


{Xi - XjY + CTj^ 


*L(x„cri) -J^Kir^Riixj) 

1=0 

- K{xi,ai,Xj), 


(19) 

( 20 ) 


where ctr and uj lie on the same grid (xi). We do not prove that iffj’^ is strictly positive definite, since the algorithm 
has to diverge if this is not the case and we have to check the convergence numerically anyway. 

The iteration is started with the initial function 




Ro{x) = —L{x,ai) 


( 21 ) 


which is motivated by the following definition of the Dirac-Delta-function 


S{x) = lim — ■ 


o-I^O TT x^ + af 


2 ■ 


( 22 ) 


Let us remind that extracting R using the limit (Ti ^ 0, i.e. 


(Jl (Jl 

R{aji) = lim —L(crR,(Ti)= lim — / dw if (ctr, cti, a>) i?(a;) . 

o-I^O TT CTi^O TT Jo 


(23) 


is numerically unstable. 

The smallest characteristic number of the integral kernel if is the inverse of the largest eigenvalue yi of {K^p'^). In 
this work we have used A = —. This value is not optimized. The larger A is, the faster the algorithm will converge, 

Ml 

but the procedure will also become more unstable. 

Under real conditions one does not know the solution of the integral equation. One has to compare the “input” 
Lorentz L with the LIT of the actual result for i?;. The relative error at every grid point 


£i{xi) 


- L{xi,ai) 

L{xi,(Ji) 


(24) 
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after I iterations and the mean square sum of ei(xi) at every grid point 


1 

(25) 

i=0 

is used to check the quality of the iterated solution Bi after I recursions. We remind here that this problem is ill-posed 
and therefore the solution is unstable under perturbation of the input data. With increasing number of iterations the 
solution will become more and more unstable. Therefore one has always to check whether the result Ri is numerically 
still stable, i.e. free of unphysical oscillations. Solving ill-posed problems requires a regularization scheme. This is 
naturally provided by this method through the grid (xi). 

Let us remark at this point that the simple integration formula (eq. 11811 1 has the advantage that all information 
entering the inversion process has the same weight except at the boundaries and it has constant grid gap, which 
provides a constant regularization over the whole integration region. In contrast, a Gauss-Legendre integration grid 
has smaller grid gaps at the boundaries than in the center of the integration interval leading to a weaker regularization 
at the boundaries and likely causing a lower quality of the inversion close to threshold. In addition the integration 
kernel K (eq. (EnH will become more easily numerically singular for the case that L is given on a grid with many 
points, i.e. more than about 50 grid points. On the other hand if L is given on a sparse grid, the Gauss-Legendre 
integration is preferred, since it provides a more precise integration for the same amount of grid points than formula 
m- For such a case the mentioned disadvantages of the Gauss-Legendre grid regarding regularization and numerical 
singularity are not so important any more. 

The iteration ca is stopped if £i is smaller than 2.5T0 if £i or maXjgjQ ^ ^„}{£;(a;i)} is increasing or unphysical 

oscillations arise in the solution. We allow at first some “free” iteration steps, depending on the quality of the input 
data, before we check if the iteration has to be terminated. This is necessary, since in the first iterations the errors £i 
or maxig{o,i,...,rt}{£/(a;i)} may increase temporarily due to large initial changes in Ri. 


D. Inversion using Banach’s fix point theorem 

Let us first consider notations and conventions used in this subsection. We are working in the Hilbert space L 2 of 
the square integrable real functions R{ijj) = {x!\R) with real arguments uj which fulfill 


J doj \u}){ui\ = 1 . 

The scalar product and the norm are defined as usual: 

/ OO nOO 

dw (i?i|w)(a;|i?2) = / dw i?i(w)i?2(w) , 

-OO J —OO 


(26) 


(27) 


and 


||i?|| := ^/(ip) < OO . (28) 

For a fixed ai > 0, we are defining the Lorentz integration operator L(ai) (LIT) in accordance with as follows: 

{u;\L{ai)\R) =L{u},ai) := f dw'- - L , > (29) 

J (UJ - uj'Y + erf 

where in general, if not noted differently, the upper and lower limit in the integrals is 00 and — 00 , respectively. For 
the inversion method proposed in this subsection, let’s introduce for a fixed Lorentzian L(w, cti) and a fixed cti > 0 
the following mapping acting on an arbitrary function Ri £ L 2 '. 


{uj\f(ai,L)\RY 


f(ai,L){RYuj)) 


—L(uj,ai) 

TT 


cri 


dw 


, Rx(uj')-Ri(uj) 
(uj' - ujY + uf 


TT 


(30) 
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As will be shown in Appendix ixn this mapping has as fixpoint just the desired response function i?(w) which 
creates via (1^ the input Lorentzian L{u!,ai): 


f{ai,L){R{uj)) = R{uj) 


R{lo) = —L{uj,ai)-— f dw'^ 

TT TT J (W 


, i?(w') - R{uj) 

- Lu)^ + af ■ 


(31) 


Now, in appendix I A 21 we will prove with the help of Banach’s Fixpoint Theorem that there exists an e > 0 so 
that for 0 < (Ti < e the series (n > 0, S L 2 arbitrary) 


has a unique fix point, which is therefore just the searched response R{u!): 


(32) 


R^^\uj) —> R^^{uj) = R{uj) for n —» 00 , 0 < cti < e . 


(33) 


This method works therefore as follows: For a given Lorentzian L{uj,(Ti) with ai sufficiently small, we choose some 
arbitrary starting function (in our practical applications, we have taken always = 0 for the sake of simplicity) 
and then calculate the series till we have reached convergence. The found fixpoint function R^^ is then identical 
with the desired response function R. This method is therefore completely parameter free and does not need any 
suitable basic functions. 

In the version discussed so far, this approach has however one serious drawback. Due to the iterative procedure, it 
is obvious that the fixpoint R^^ is a function, i.e. differentiable up to infinite order. In other words, this method 
is strictly valid only for response functions which are (7°°. The physical response function, however, is not analytic 
at the threshold energy wth of the considered inclusive reaction and is therefore not C°°. In order to get rid of this 
problem, one has to modify slightly the above discussed algorithm. Instead of the range ] — 00 , (X)[, all occurring 
functions L, R,... are now only defined on the interval [a, 6], where a = c^th- The upper value b should be in principle 
00 . However, in practical applications, the Lorentzian L, which serves as a input in the algorithm, is known only in 
some finite range of w, so that in practice b is not infinite, but still “large”. 

Instead of 03) we have therefore to consider 



and the counterpart of the linear mapping for finite integrals becomes then 


with the function 


The series 


T(a,b){cri, L){Ri (w)) 


cri 


0‘{a,b) i^) 


CTl 


^(a,b) (WjCn) 


^{a,b) (^) 




(w' - w)2 + af 


, , oj — a. ,b — uj 

a(a,b)i^) = arctan(-) + arctan(-) 




cri 


R 


(n+1) 

(a,fc) 


(w) :—f(a,b){'yi,Lfa,b)){R, 


\a,b) 


W)) 


(35) 


(36) 


(37) 


has again, as will be shown in a.nnendix IA 31 the desired response function R as unique fixpoint for sufficiently small 
Cti. So the method works also for functions which are defined only on a finite range. 

Summarizing, we have proven that this method allows for not too large cti a unique determination of the response 
function R for a given Lorentzian L. As discussed in section nvi in detail, it may however happen in practice that L 
is not very precisely known, containing for example oscillations due to numerical uncertainties or approximations in 
its explicit calculation. Therefore, similar to the Fridman approach, also in this method a sort of regularization has 
to be used in order to avoid resulting unphysical oscillations in R. In practice, two numerical tools are used which 
mutually supplement each other, leading both to a reduction of the numerical resolution. First of all, the number 
of mesh points in the occurring integrals can be reduced (in analogy to the Fridman method) and second, the total 
number of iterations should be restricted, for example by choosing a suitable termination condition for £i (ESI . 









IV. DISCUSSION OF RESULTS 


In the following we apply the various inversion methods of section III to two different cases. First we investigate 
the efficiency of the regularization procedures implemented in the various inversion methods. As test case we consider 
a realistic example, where L((Tr, cti) is calculated from eqs. 1111511 without any knowledge of R{uj). To this end we take 
L{a^, ai) from ref. where the longitudinal response function for inclusive electron scattering off the deuteron at a 
momentum transfer of = 5 fm“^ is considered. As described in ref. eq. © was solved in an approximate way 
thus leading to an L{aB.,cri) with a non-negligible error margin. It is evident that for such a case a regularization 
must be present in the inversion method otherwise reasonable results cannot be obtained. Our second example serves 
for a different test, namely to check the ability of the various inversion methods to precisely reproduce response 
functions with rather complicated structures. For this high precision check we take a test case, where the response 
function R{uj) is known beforehand analytically, and choose a function with a double peak structure. Then T(ctr, cti) 
is calculated via eq. CJ with a rather high numerical precision and, finally, the obtained LIT is inverted. 

The LIT of our first case is illustrated in fig. |2 One sees a pronounced peak at about 50 MeV and one notes slight 
oscillations of the transform at higher energies. As already mentioned the LIT is taken from P|, where the longitudinal 
response function (?) for inclusive electron scattering off the deuteron at a momentum transfer of = 5 fm“^ 

is considered. In ref. [l| eq. © was solved with intention only approximately to obtain a transform with some error 
margin in order to work with a test case of not too high numerical precision. For example the visible high-energy 
oscillations of the LIT in fig. |21 are due to computational approximations. The solution of © was obtained in [y 
by putting 'l'(r) equal to zero at a neutron-proton distance of r = 20 fm instead of using a more exact asymptotic 
behaviour. Of course in case of the deuteron a more precise LIT could be calculated easily. 

We would like to emphasize once again that in case of such a not very precisely determined LIT a proper regular¬ 
ization must be built into the inversion method otherwise reasonable results cannot be obtained. In [J the standard 
inversion method was used for the inversion and problems were encountered, since the results exhibited oscillations 
at lower energies. The problems arise from Coulomb monopole (CO) and quadrupole (C2) transitions, which lead to a 
low-energy shoulder in the response. As will be discussed also later on the standard inversion method has difficulties 
to reproduce a more structured response function. In P| these problems could be overcome by making separate 
inversions for the CO and C2 LITs. 

First we discuss the inversion results with the wavelet method. Here we should mention that Rl contains also an 
elastic contribution. Thus it is necessary to introduce an additional 5-shape basis functions which accounts for it. 
This can be achieved easily and the inversion result leads not only to the inelastic response function, but also to a 
value of the elastic form factor at = 5 fm“^. In principle one could also calculate the elastic form factor separately 
and then subtract its contribution from the LIT. In fig. 0 we show two inversion results. While the inversion with 
the parameters M = 0 and A’q = 15 leads to a smooth curve, one already finds strong unrealistic oscillations if Nq is 
increased to 31. Even stronger oscillations are found for the inversions with the other parameter values of table 1. It 
shows that for the present example a rather low resolution, i.e. a rather strong regularization, is necessary in order 
to obtain an Rl result which can be interpreted as realistic. 

In contrast to the wavelet method the Banach inversion and the Fridman method are iterative methods. The nature 
of the iterative inversion method is that with increasing number of iterations the error £i lea. (l25ll 'l will become smaller 
and smaller but the solution will become more and more unstable. That means that the amplitudes of unphysical 
oscillations will increase. In the following the integration grid is used as regularization. Therefore the larger the grid 
gaps are the smaller the frequencies of the unphysical oscillations will be. In principle one could also introduce a 
regularization operator and calculate the optimal number of iterations depending on the desired error £i lea . 12511 1. 
We did not follow this approach of regularization since the grid-regularization works successfully as discussed in the 
following. To simplify the inversion process the elastic contribution to the LIT, Lei = 0.0436/((£’o — ctr)^ + c^?)) with 
—Eq = Lb = 2.2246 MeV and cti = 10 MeV Q, has been removed in both inversion methods. 

In fig. 0 we illustrate examplarily the procedure to determine the inversion for the Fridman case. First the number 
of iterations is fixed to 120 and the integration is done on a grid with 40 points (curve: pts:40, it:120). This curve 
shows still quite large unphysical oscillations although the integration grid consists of only 40 grid-points. By further 
reducing the amount of iterations to 20 (curve: pts:40, it: 20) the amplitude of the unphysical oscillations is reduced. 
Then we lower the amount of grid points till a smooth curve is obtained. At that point a strong regularization is 
introduced. To improve the threshold region the grid (xi) is split up into two grids, a small inner grid using the 
simple integration method (eq. (I18II 1 and an outer grid using a Gauss-Legendre integration grid. Finally a systematic 
variation of the number of grid points, grid boundaries and the number of iterations yields a parameter set where 
the result is stable under small changes of these parameters. The result shown fig. 0 and fig. 0 is obtained with 12 
iteration where 6 points are used for the inner grid [0, 30] MeV and 14 points for the outer one ]30,197] MeV. 

A similar pattern occurs also for the Banach method. Concerning the chosen resolution, we use here 24 Gaussian 
mesh points for the evaluation of the integrals and stop the iteration process for £i < 10“®. In contrast, if the 
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Lorentzian is numerically better determined, like in our second example discussed below, typically a couple of hundreds 
of mesh points and termination conditions like £i < 10“^^ (or even less) can be used for small aj without producing 
unphysical numerical oscillations in the resulting response R. 

Our deuteron example has the great advantage that -R(oj) can also be determined directly from eq. Q, since two- 
nucleon continuum wave functions can be calculated without problems. In fig. I^we compare this response function, 
which is also taken from Q, to those from the three inversion methods. It is evident that in spite of the strong 
regularization the inversion results lead to an overall correct description and form a kind of error band. The wavelet 
and Banach inversions oscillate around the true response with an error of less than 2 or 3 % in the peak region. 
The Fridman inversion exhibits a different pattern, but also describes the peak region quite well. The relative errors 
become somewhat larger at lower and higher energies. Nonetheless in comparison to the standard inversion method, 
which leads to very strong low-energy oscillations (see Q) , our present results have a much smoother and more correct 
form at lower energies. It is important to note that errors of the inversion results are due to a not sufficiently precisely 
determined LIT. In fact as shown in ref. Q the error of the LIT amounts up to 1.5 % in the peak region and up to 
about 4 % beyond 120 MeV. One sees that an error in the transform leads to an error of similar size in the inversion. 
Thus more precise results could be obtained by calculating the LIT with a higher precision. However, even with the 
present precision the low-energy shoulder could be described much better if one makes separate inversions for the CO 
and C2 transition strength as was done in ref. P|. 

After our first example, which checked the proper implementation of the regularization into the various inversion 
methods, we turn to our high-precision test. In fig. Owe show the chosen i?(a;) up to w = 120 MeV. The response 
function exhibits two peaks, one at w = 20 MeV and the second at w = 50 MeV. Both peaks have the same heights, 
but different widths. The corresponding LITs with aj = 5, 10, 20, and 40 MeV are illustrated in fig. [3 It is evident 
that only the smallest value of ai leads also to two peaks for the transform, while the information about the two 
separate peaks is more and more smeared out with increasing crj. 

In order to invert the LIT we calculate L(crR, (Ji) for each of the four cri values from ctr = —10 MeV to ctr = 210 
MeV at 441 equidistant grid points. It is instructive to use the standard LIT inversion method, described in section 
IIIA, for the transform with the highest resolution (tri = 5 MeV). In fig. I^it is shown that the result is not of a 
very high quality. There are low-energy oscillations and the two peaks are not well reproduced. One sees that the 
standard inversion method has difficulties to reproduce a more structured response function like our present double 
peak example. The reason for the problem are the basis functions. They have a too large extension in the w-space, 
therefore oscillations are introduced easily. In case of a more simpler structure like a single peak response such 
difficulties do not appear (see e.g. El)- We do not show results with the standard inversion method for larger values 
of ai- As one may expect they lead to even worse results. 

It is obvious that alternative inversion methods are needed in case of more structured response functions. In fact 
the three inversion techniques described in sections IIIB-D lead to much better results. For ai = 5 and 10 MeV one 
obtains results so close to the true R that in fig.|Slthey would all appear identical to the true response function. Thus 
in fig. 13 we show the relative errors of R in an w-range, where R{oj) is not too close to 0. It is evident that the relative 
errors are very small: generally much less than 0.001 for ui = 5 MeV and less than 0.01 for cti = 10 MeV, only close 
to the low- and high-energy borders, where R approaches 0 the error can become a bit larger. Here we would like to 
point out that the threshold region can also be inverted with greater care (see final paragraph of subsection IIIA). In 
figs, uni and CH we illustrate the inversion results for the two higher ui values. One sees that one still obtains rather 
good results with ui = 20 MeV, particularly with wavelet and Banach’s fix point theorem inversion techniques. With 
CTi = 40 MeV the inversion results become worse. The shapes of the two peaks are not very well reproduced, but at 
least the peaks are recognized as two separate structures and also the high-energy tail is still rather well described. 
However, it is obvious that for the present case a ai value of 40 MeV is not sufficient to have precise inversion results. 
On the other hand with a further increase of the numerical precisions in calculating the transform and the inversions, 
which is nothing else than a further reduction of the regularization, one could have also for ai = 40 MeV much better 
results. 

In the above example Fridman and Banach methods lead to a single inversion result for a given cti, while the wavelet 
technique leads for any set of values for M and Nm (see table I) to principally different inversion results, in absence 
of unphysical oscillations they are of course almost identical. For ai = 5, 10, and 20 MeV the result with the best 
description of L(tTR,CTi) (sum of quadratic errors) is taken (cti = 5 and 10 MeV: M = 2 and Nq = Ni = N 2 = 63, 
Cti = 20 MeV: M = 2 and Nq = 63, Ni = 31, N 2 = 15). For ui = 40 MeV, however, the inversion result with the 
smallest error shows strong and unrealistic oscillations (see fig. I12II and is thus discarded. Also the 14 next best fits 
lead to similar unrealistic oscillations showing that also for this case a rather strong regularization is necessary. Only 
the 16th best result (M = 0, Nq = 63) does not exhibit such an unrealistic pattern and is therefore shown in fig. 1111 
Such strongly oscillating solutions like that of fig. El can be identified easily as unrealistic, since a true response 
function has to be positive definite and in addition it is difficult to imagine that a true response can exhibit such a 
regular oscillation pattern. 
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We summarize our results as follows. The various newly introduced inversion techniques for the Lorentz Integral 
Transform contain a proper regularization method and are thus able to treat the ill-posed inversion problem. In 
addition all these techniques are capable to invert also transforms of responses with rather complicated structures 
with a very high precision. For such cases they lead to considerably better inversion results than the standard LIT 
inversion method. Our results show that the LIT approach is not restricted to treat cases, where simple structures 
appear, but is suitable for many different applications. 


APPENDIX A: MATHEMATICAL PROOFS CONCERNING THE PINPOINT METHOD 

1. Proof of (1311) 

In order to prove iU, let’s consider the expression 

, Ric') 


B(Lo,ai) := 3 / dw' — 

J w' 


— u! — lai 

which can be rewritten in two different manners: First, one has the identity 

, R{lo'){uj' — uj + iai) 


B{lo, ai) = SJ y du! 


and second, one can write 


- Lu)^ + af 

,R{uj') - R{uj) 


= CTiL(w,cri) 


B{io,ai)=Q [ duj'^^^ — ?^ + ^(r{uj) [ dJ - ) 

J u' — UJ — l(7i \ J Uj' — UJ — IG\ ) 


— UJ — IGi \ J UJ' — UJ — IG\ ^ 

The second integral on the right hand part of (IA3II can be easily solved with the help of 

1 


dx 


+ a? 


= 1 , (a > 0) . 


(Al) 


(A2) 


(A3) 


(A4) 


Comparison of (IXali and (EH) yields then m- 


2. Proof of uniqueness of iteration procedure 

In this appendix we prove that there exists an e > 0 so that for all ai < e the series (EU converges uniquely against 
a function R^^ which is identical with the desired response R{ijj). 

For that purpose, let us introduce the difference 

:= i? - R^^'>. (A5) 

Exploiting P|l. ljiTT|l .(P)l. the function A fulfills 

A(’"+i) := TA(") (A6) 


with the linear mapping 


TT J [uj' — UjY -I- af 


(A7) 


At next, we will prove that T is - for sufficiently small (Ti - a contractive linear mapping, i.e. there exists a constant 
0 < g < 1 so that 


||TA(")|p = ||A('*+i)|p < g||A(^ 


■)||2 


(AS) 


Proof: 
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For the norm of we obtain 

||A(”+i)|p = J da;A("+i)(u;)A(”+i)(a;) 


= —^ duj dio' I duj" 


A(") (w') - A(’^) (w) \ / A(”) (w") - A(”) (w) 


{lo' — loY + ct^ ) \ {lo" — ujY + cr^ 


(A9) 


which can be rewritten as 


iiA‘”+'>it = ^ J J J 

(A('‘)(a;')A(")(w") + A(")(a;)A('‘)(w) - A^"’“ A(")(a;")A(")(a;)) . 


- loY + arf 


(uj" - u;Y + 


(AlO) 

This expression can be simplified with the help of (IA4II and with the identity (which can be easily proven using the 
method of residues) 


dhJ 


(w' — tt’)2 + (Tj j \{ijj" — ujY + J O’! \ (w'— + 4crj 


27r 


as follows: 


||A("+i)||2 ^ ||a(”)|| 2 + ^(A(")|L(2fTi)|A(")) - — (A^|L(ai)|A('‘)) . 


For further exploitation, let’s prove at next the following 
Lemma: There exists an e > 0, so that for all cti < e 

2(Ti 


— (A(’")|L(2cri)|A(”)) < A(cri)(A(”)|A(’")), 


where the upper limit A(cri) for ctj < e fulfills the estimate 


1 — 5<A<l + 5, 5 < - . 

O 


Proof of the Lemma: 

Due to 


CTi 


lim „ 

o-i^oo TV [x — uY + CTf 


= 5{x - y) , 


the statement is obviously true for cti —> 0 with 


lim A[(7i) = 1 . 

(Ti—^0 


(All) 

(A12) 

(A13) 

(A14) 

(A15) 

(A16) 


Because the LIT-operator in (1^ is depending analytically on the parameter cti, the statement follows immediately. 
With the help of the Lemma, one obtains for the following estimate: 


||A('*+i)||2 ^ |||A(")||2 + ^(A(”)|L(2cti)|A(’^)) - —(A(”)|L(cti)|A(’^))| 

TT TT 

< |l + (l + ,5)-2(l-5)| ||A(")||2 

= 3^||A(”)||2 . (A17) 

Due to 

g := 3(5 < 3^ = 1 , (A18) 

O 

the mapping A^”) ^ A^”+^) defined via (lAbll is therefore for ai < e contracting. 

In consequence, we can now apply Banach’s Fixpoint Theorem: It states that the series A^"^ converges uniquely 
against a fixpoint A. Obviously, (IA6II is fulfilled by A^"^ = 0 for all n, so that therefore 

A(") ^ A = 0 (A19) 

and therefore, due to 1A5II . ^ R uniquely. 
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3. Incorporation of correct threshold behaviour of the response function 

If the lower and upper value of the integral ll^ are not any longer ±oo, but some finite values a and b, the proofs 
presented in aopendix IA ll and IA 2l have to be solely repeated where instead of the expression (IA4II its counterpart for 
finite integrals 


^(a,b) (^) 


dx 


1 


a;2 + 


= 1 


(a > 0) 


has to be used (the function a is given in 13(111 1. Therefore, instead of one has the identity 


(A20) 


R{u}) = - ^—-—L{uj,(7i)— ^ 


0:{a,b) 


0^{a,b) (w) 




(w' - w)2 + 


2 ’ 


(A21) 


which motivates therefore the mapping L) according to 13511 . 

The counterpart of the difference function 1A5II now becomes 








a(a,b)(w) 


duj 


/ (a 




( lu ' - cc)2 + 


The norm of this function 


J a 


= dco duj' 


dio" 


O’! 


(a,6)(w)y +^1 / \iuj” - Uj)"^ + 






(A22) 


(A23) 


is much more complicated to be handled than its counterpart ITinll : Due to the finite limits in the integrals, the 
method of residues cannot any longer be used. Moreover, the constant factor tt in (rnm now turns into the complicated 
function Q!(a,h) (w) so that a simple relation like 1A12II is not any longer valid. However, recalling the representation 
(Tn^ of the j-function as well as (for arbitrary a, 6; provided a < uj < b) 


lini a(a,b)(w) = tt 


(A24) 


it is obvious that for sufficiently small cti the influence of the finite limits a and b in the integrals dies out so that 


again the mapping is a contractive one. 


,(«■+!) 
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FIG. 1: Mexican hat wavelet (eq. (10)). 



FIG. 2: LIT of the first test case taken from Q 







R [MeV 1] •• R [MeV-1] 



FIG. 4: Fridman inversion results for the first test case (full: 20 grid-points and 12 iterations, dashed: 40 and 20, dotted: 40 
and 120). 
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FIG. 5: Final inversion results with the three new methods for the first test case are shown (upper panel; full: original response 
function, dashed: wavelet, dashed-dotted: Fridman, dotted: Banach) along with their relative error (lower panel). 



FIG. 6: Double peaked response function of second test case. 
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FIG. 7: LIT of the double peaked response function with different ai (curves for cti = 10, 20 and 40 MeV are scaled by factors 
of 2, 4 and 8, respectively). 



FIG. 8: Result for cri = 5 MeV (dashed) obtained with the standard method (sec. IIIA) for the second test case and the 
original double peaked response function (full). 








FIG. 9: Relative error of the inversion results for the second test case obtained with the new methods (full: wavelet, crosses: 
Fridman, points: Banach) for ai = 5 MeV (upper panel) and for ai = 10 MeV (lower panel). 








18 




FIG. 10: Results of the three new inversion methods (full: original response function, dashed: wavelet, dashed-dotted: Fridman, 
dotted: Banach) for the second test case with cti = 20 MeV are shown (upper panel) along with their relative error (lower 
panel; full: wavelet, crosses: Fridman scaled by a factor of 0.2, points: Banach). 
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FIG. 11: Results of the three new inversion methods (full: original response function, dashed: wavelet, dashed-dotted: Fridman, 
dotted: Banach) for the second test case with cti = 40 MeV are shown (upper panel) along with their relative error (lower 
panel). 















